Anomalous Spreading of Power-Law Quantum Wave Packets 



o 

CD 

Q 

(N 



Fabrizio Lillo and Rosario N. Mantegna 
Istituto Nazionale per la Fisica della Materia, Unita di Palermo 

and 

Dipartimento di Energetica ed Applicazioni di Fisica, Universitd di Palermo, Viale delle Scienze, 1-90128, Palermo, Italia 

We introduce power-law tail quantum wave packets. We show that they can be seen as eigenfunctions 
of a Hamiltonian with a physical potential. We prove that the free evolution of these packets presents 
an asymptotic decay of the maximum of the wave packets which is anomalous for an interval of the 
characterizing power- law exponent. We also prove that the number of finite moments of the wave 
packets is a conserved quantity during the evolution of the wave packet in the free space. 
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Power-law probability density functions W are receiv- 
ing a lot of attention in different research fields ^l). 
Stochastic processes with power-law distributions may 
or may not be characterized by a typical scale in time 
and/or in the size of the random variable. One implica- 
tion of the absence of a typical scale is the divergence of 
the variance of the distribution. Examples of phenomena 
without a typical scale are observed in physical systems 
at the critical state in self-organized || and in com- 
plex systems || . Considering few-body quantum systems 
- (i) scale free power-law processes have been observed in 
a quantum system by investigating and modeling exper- 
iments of velocity selective coherent population trapping 
@,| and (ii) the power-law temporal growth of the mo- 
ments of a wave packet has been investigated in quantum 
systems with fractal energy spectra and eigenfunctions 
such as the Harper model ■ Spatial power-law wave 
functions have not been considered within the framework 
of non-relativistic quantum mechanics. The probabilis- 
tic interpretation of the wave function and the recent re- 
sults on power-law distributions observed in physical sys- 
tems motivates us to investigate the properties of power- 
law wave functions in quantum mechanics. In this letter 
we consider quantum wave packets with power-law tails. 
Specifically we focus on - (i) their physical properties. 
Namely the uncertainty product, the associated energy 
and the momentum distribution and (ii) the spreading of 
such quantum wave packets during the free evolution. 

We define as Power-Law Tail Wave Packet (PLTWP) a 
wave function ip(x) describing a non-relativistic spinless 
particle in one dimension, which decreases with x as 



tp(x) 



(1) 



This class of wave packets is square-integrable only if 
a > 1/2. For the sake of simplicity, in this study we as- 
sume that wave function is real, positive and even. The 
study of quantum wave packets with wave function real 
or complex, uneven and with zeros is presented elsewhere 
pT| . One of the properties of power-law distributions is 
that only a finite number of moments of the variable are 
finite. In quantum theory this implies that the moments 



of position operator (x m ) with m > 2a — 1 are infinite. 
The lack of finite moments of x has a counterpart in the 
properties in k — of Fourier transform g(k) = FT[ip(x)] 
which gives the amplitude probability distribution of mo- 
mentum. We note that when n < a < n + 1, only the 
first n — 1 derivatives of g(k) exist in k = 0. In the special 
case 1/2 < a < 1 the Fourier transform g(k) is infinite in 
k = 0. The behavior of g(k) when k ~ is described by 
two different series expansions depending on the value of 
a. Specifically, when a is not an odd integer number 

g(k) ~ a + a 2 k 2 + ... + a 2n k 2n + b\k ^ +o(| k | Q ), 

(2) 

where 2n is the largest even integer number smaller than 
a — 1 . When a — In + 1 the series expansion is 

g(k) ~a + a 2 k 2 + ... + a 2n k 2n + bk 2n log | k | +... . (3) 

The case a = 1 has a logarithmic divergence in k = 0. 

In spite of this behavior we wish to point out that 
all moments of the momentum operator (including the 
kinetic energy) of the particle (p m ) (m > 1) are finite. In 
fact, a property of Fourier transform is 



FT[^ m \x)] = {-ik) m g(k), 



(4) 



where ip^ m \x) indicates the m-th derivative of ip(x). 
This property is valid under the hypothesis that 
limui-foo ip( r \x) — for r = 0, l,...,m— 1. Since the 
tails behave as | x | _Q ~ m ; -0( m )(x) is absolutely intc- 
grable. By using Riemann-Lebesgue lemma p"2[ , one can 
conclude that lim^oo k m g(k) = for all m. 

A consequence of the finitcness of momentum moments 
is that the uncertainty in momentum Ak is always finite 
for PLTWP. As mentioned above, when a < | the sec- 
ond moment of the position operator and thus the root 
mean square deviation Ax of position are infinite. There- 
fore this kind of wave packets have an uncertainty prod- 
uct AxAk which is infinite when a < |. This physical 
property reflects the fact that for PLTWP with a < | a 
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typical scale exists in momentum space whereas the wave 
packet is scale free in space. 

Can a PLTWP be an eigenfunction of an Hamiltonian? 
To answer this question we note that each PLTWP can 
be related to a specific physical potential. The PLTWP is 
then an eigenfunction of the corresponding Hamiltonian. 
The equation providing the potential associated with a 
particle described by a wave function ipix) is |l3|| 



U(x) = E 



1 d 2 ip(x) 



2M ip(x) dx 2 



(5) 



where E is the eigenvalue of energy and M is the mass 
of the particle. The shape of the potential depends on 
the local properties of the wave function. ^From Eqs. 
(1) and (5) it is immediate to conclude that the poten- 
tial associated with PLTWP behaves asymptotically as 
x~ 2 [Q. As an illustrative example we present PLTWPs 
defined as 



il>{x) 



N 



{x 2 +7 2 )"/ 2 ' 



(6) 



where TV is a suitable normalization constant and 7 is 
a scale parameter. The above family of quantum wave 
packets is related to the Student's ^-distribution when a 
is integer. The associated potential is 



_/^_ x 2 (l + q)- 7 2 
[X) ~2M a (x 2 + 1 2 ) 2 ' 



(7) 



and is shown in Fig. 1. In this figure the eigenvalue 
E is set equal to zero. The potential is a single well 
potential with two symmetrical confining barriers. The 
potential reaches a maximal value and then decreases 



asymptotically as U (x) 



A general property of this 



potential is that by increasing the value of a the depth 
of the potential well increases. From Fig. 1 it is clear 
that the associated potential does not present anomalies 
of any sort. 

Hereafter we consider the properties of a PLTWP in 
the simplest case of dynamical evolution, namely the free 
wave evolution. We assume that a t — the wave packet 
in the free space has the asymptotic properties of Eq. 
(1). We focus our attention on the spreading of the wave 
packet. During the dynamical evolution in free space the 
wave function at time t is given by 
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g(k)e^ kx - Mk2 / 2M Uk. 



(8) 



We briefly recall that in the case of a Gaussian wave 
packet the amount of spreading of the wave packet can be 
quantified cither by considering the time dependence of 
the position variance or by determining the time depen- 
dence of the maximum of the wave function. To quantify 
the amount of spreading of a free wave packet in the same 
way for PLTWP with finite or infinite variance we choose 
to focus our attention to the asymptotic behavior in time 



of the wave function in a specific position (for example 
x = 0). In the Gaussian case, the variance of | tp(x,t) \ 2 
is asymptotically proportional to t 2 and the maximum of 
the wave function | ip(0,t) \ 2 decreases as t -1 asymptoti- 
cally. We will show in the following that this asymptotic 
behavior observed for Gaussian quantum packets is not 
universally observed in the free evolution of a quantum 
wave packet. 

We consider here the free evolution of the maximum of 
the packet, which is described by Eq. (8) with x = 0. It 
is possible to give an asymptotic expansion in time t of 
integral of Eq. (8) using phase stationary method (see for 
example |jq]). This method shows that the asymptotic 
expansion of ip(0, t) is determined by the dominant term 
of series expansion about the origin (k = 0) of g{k) and 
by the function —htk 2 /2M. The first term of the series 
expansion about k = of g(k) is g(k) ~ Q|fc| A_1 for any 
value of a ^ 1. By using phase stationary method [fTsf , 
we know that when < A < 2 it is possible to give an 
asymptotic expansion of integral of Eq. (8) with x = 
in the form 
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2tt 



(9) 



where (3 = ht/2M . If g(k) has a finite non- vanishing limit 
as k goes to zero then A = 1 and we conclude that the 
maximum of the packet decreases asymptotically in time 
as l/yt. This is the case, for example, of Gaussian wave 
packets cited above and, in general, of all commonest 
wave packets with finite FT in k = 0. From Eqs. (2,3) 
and (9) we conclude that this is also the case of PLTWP 
with a > 1. These packets show a customary asymptotic 
dynamics of the maximum. PLTWPs with 1/2 < a < 1 
show a different behavior. In fact, in these cases g(k) 
diverges in k — as | k \ a ~ 1 and, as a consequence 
A = a < 1. Hence, form Eq. (9) we conclude that 
the maximum decreases asymptotically as t~ a l 2 . This 
is a maximum decrease, which is anomalous and slower 
with respect to the decrease observed in customary wave 
packets. The i -1 / 2 behavior observed in Gaussian wave 
packets is interpreted in terms of the time evolution of a 
group of classical particles with a momentum dispersion 
Ap. A similar simple picture cannot explain the behavior 
observed for PLTWPs with a < 1 @]. 

To illustrate the process of convergence of different 
wave packets to the expected asymptotic behavior we cal- 
culate the time evolution of the amplitude of wave packet 
at x = for three different cases. Specifically we consider 
a Gaussian wave packet and two PLTWPs of the class de- 
fined by Eq. (6) with a = 3 and a = 0.75. In Fig. 2 we 
show the numerical estimate of \ip(0, t)\. In the figure is 
clear that the Gaussian and the PLTWP with a = 3 
soon converge to the usual 1/Vt asymptotic behavior 
whereas the PLTWP with a = 0.75 slowly converges to 
the anomalous asymptotic behavior of l/W 2 = 1/i 375 . 

The case a = 1 cannot be handled with phase station- 
ary method because of logarithmic divergence of g{k) in 
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k = 0. Although we are not able to provide a general 
answer for the case a = 1, we are able to determine the 
asymptotic behavior in the specific case of a packet de- 
scribed by Eq. (6) with a = 1. For a such wave packet, 
the maximum decreases as 1/ y/i. 

In order to obtain a more complete description of the 
free wave spreading, we consider the time evolution of 
the tails of the packet. We prove that the free wave 
packet evolution of a PLTWP conserves the tails. More 
precisely, if at t = the dominant term of asymptotic 
expansion of tp{x) is c | a; |~ Q , at each subsequent time 
the asymptotic expansion of ip(x, t) will be dominated 
by c | x \~ a . In this sense, the free evolution cannot 
change the asymptotic properties of the packet. We can 
visualize this property dividing the whole set of PLTWPs 
in classes such as each class is characterized by a value 
of a. In other words, any packet belongs to the same 
class at any time during the free wave evolution. Since a 
determines the number of finite moments of operator x, 
our result implies that position moments cannot became 
finite if they were infinite at t = and vice versa. We 
prove our statement as follows. Let us first consider Eq. 
(8) and look for the asymptotic expansion in x of the 
FT of f(k) — g(k)e~ % @ k . The asymptotic expansion of 
the FT of f(k) is determined by its singularities Jl2| . 
By following Ref. fl2]| , we say that f(k) is singular in 
fco if one cannot differentiate f{k) in ko any number of 
times. From properties of g(k) of PLTWPs we conclude 
that k = is the only singularity for /(fc). In order to 
find the asymptotic expansion of Eq. (8), we construct 
a function F(k) in such a way that f(k) = f(k) — F(k) 
has absolutely integrable m-th derivative in an interval 
including k — 0. F(k) must be a linear combination of 
powers of k and product of powers of k and logfc [ fl2| . 
Moreover the m-th derivative of f(k) must be absolutely 
integrable in an interval from a real value to infinity. If 
all these hypothesis are verified, the FT of f(k) is equal 
to the FT of F(k) plus o(| x |~ m ). Depending on the 
specific value of a the regularizing function assumes a 
different form. In spite of this, our results do not depend 
on the specific value of a. For the sake of simplicity, we 
present here the case with In + 1 < a < 2n + 2. The 
demonstration of cases for other value of a (specifically 
In < a < 2n + 1, a = 2n and a = 2n + 1) are similar. 
In order to regularize the 2n + 1-th derivative of f(k) the 
natural choice is 

/(*) = g(k)e- il3k2 - b | k r 1 . (10) 

In fact the first 2n + 3 derivatives of f(k) are absolutely 
integrable in an interval including k = 0. Therefore, we 
obtain the asymptotic expansion of ip(x, t) as 

^(x,t) = FT(F(k)) + o(\ x \- 2n - 3 ) = 

c\x\~ a +o(| x r 2 - 3 ) (11) 

which demonstrates our assertion. 

As an illustrative example let us consider the free evo- 
lution of the PLTWP of Eq. (6) with a = 2 and 7 = 1. 



The corresponding wave function has the form of the 
Cauchy distribution and it is possible to find the ip(x,t) 
analytically. In Fig. 3 we show the square modulus of 
the positive tail of the packet versus x at different times 
in a log- log plot. The figure shows that the dominant 
term of the asymptotic expansion, i.e. x~ 4 , is the same 
at any time. The figure also shows that new terms of the 
asymptotic expansion become relevant at longer times 
and the region of asymptotic convergence moves towards 
larger values of x. 

A connection between stochastic processes and quan- 
tum processes has been considered within the framework 
of stochastic mechanics |17j |. Here we note that the free 
evolution of a quantum wave packet is closely related 
to a superdiffusive (ballistic) stochastic process. When 
a superdiffusive behavior is observed in classical [l8]-pc|] 
and quantum |^lf processes, the central limit theorem 
does not apply and a general theoretical description is 
lacking. In our study we obtain two general conclusions 
in the problem of the free evolution of a PLTWP. The 
first concerns the temporal evolution of the maximum of 
the wave packet (which corresponds to the probability 
of return to the origin in stochastic processes). In the 
well-known case of a Gaussian wave packet (which has a 
fractional Brownian motion with exponent h = 1 as the 
approximately corresponding stochastic process |p2| ) the 
variance is quadratic in time and the probability of return 
to the origin is inversely proportional to the time. Sim- 
ilarly, we find that the same behavior is asymptotically 
observed for the evolution of an even and real PLTWP 
when the exponent a is greater than one. Conversely, for 
values of a within the interval 1/2 < a < 1 we observe 
an anomalous behavior of the time evolution of the max- 
imum of the wave packet. We are not able to interpret 
this unexpected behavior on a semiclassical basis. The 
second conclusion concerns the conservation during the 
quantum time evolution of the number of finite/infinite 
moments of the t — distribution. This behavior is pecu- 
liar to this quantum dynamics and could not be observed, 
for example, in stochastic processes obeying the central 
limit theorem. 
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FIG. 1. Potential profile, Eq. (7), associated with the class 
of power- law tail wave packets of Eq. (6). Here we show a 
region of a values ranging from the minimal allowed value 
a = 1/2 (excluded) to the value a = 4. In the entire interval 
a > 1/2, the potential does not show any pathological behav- 
ior. In our calculation we set 7=1. The potential U(x) is 
given in units of h 2 /2M. The PLTWP is the eigenfunction of 
the corresponding Hamiltonian with eigenvalue E = 0. 




FIG. 2. Log-log plot of the time evolution of the wave 
packet amplitude at x = for three different wave packets. 
The usual asymptotic behavior t~ L / 2 is observed for the Gaus- 
sian (diamond) and for the PLTWP defined by Eq. (6) with 
a — 3 (square). The anomalous decay predicted for PLTWP 
with 1/2 < a < 1 is observed for the PLTWP defined by Eq. 
(6) with a = 0.75 (circle). To estimate the degree of con- 
vergence, we also show the asymptotic behaviors predicted in 
the two cases as straight lines. The dashed line indicates the 
t^ 1 / 2 behavior and the solid line shows the t -0 - 375 behavior. 
In our calculation we set 7 = 1, h/2M = 1 and we chose the 
value of Ax of the Gaussian wave packet equal to the same 
quantity of the PLTWP with a = 3 at t = 0. 
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FIG. 3. Log-log plot of the positive tail of \tp(x, t)\ 2 for 
the PLTWP of Eq. (6) with a = 2 as a function of x for 
different values of the time. The bottom line refers to the 
case t = 0. From bottom to top the other lines describe 
the cases t — 10 4 , t — 10 5 , and t — 10 6 respectively. The 
dominance of the leading term of the asymptotic expansion 
at t — is also observed for longer times at larger values of 
x. In our calculation we set 7 = 1 and h/2M = 1. 
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